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We study the distribution of threshold forces at the depinning transition for an elastic system of finite size, 
driven by an external force in a disordered medium at zero temperature. Using the functional renormalization 
group (FRG) technique, we compute the distribution of pinning forces in the quasi-static limit. This distribution 
is universal up to two parameters, the average critical force, and its width. We discuss possible definitions 
for threshold forces in finite-size samples. We show how our results compare to the distribution of the latter 
computed recently within a numerical simulation of the so-called critical configuration. 



I. INTRODUCTION 



The dynamics of elastic objects driven by an external force 
in disordered media has attracted considerable theoretical and 
experimental interest during the last years QUIlt. The rea- 
son is twofold. On one hand, elastic objects in disordered me- 
dia exhibit the rich behavior of glassy systems and thus their 
study can help us to understand the physics of more complex 
systems such as spin glasses H or random field systems H. 
On the other hand, the motion of elastic objects in disordered 
media is an adequate description of many experimental sys- 
tems. One can divide these systems into two classes. The first 
class comprises periodic systems, the most prominent exam- 
ples being charge density waves (CDW) in solids. These start 
sliding when the applied electric field becomes large enough 
01 • Vortex lines in disordered superconductors form a quasi- 
ordered periodic Brag glass phase Ola]. The second class 
includes propagating interfaces such as domain walls in mag- 
netically or structurally ordered systems KJK interfaces be- 
tween inmiscible fluids in porous media flOll or dislocation 
lines in metals |11]. To unify the mathematical description 
of these different systems one uses the notion of "manifolds". 
In all these systems the interplay between quenched disorder 
and elasticity leads to a complicated response of the system 
to an applied external force. At zero temperature, a driving 
force / exceeding a certain threshold value f c is required to 
set the elastic manifold into motion. This depinning transition 
shares many features with critical phenomena 1 12]: character- 
istic lengths diverge close to the transition as £ ~ (/ — f c )~ u 
and the system becomes extremely sensitive to small pertur- 
bations. Following the description of standard critical phe- 
nomena, one can identify the ordered phase with the moving 
phase with force / > f c , and the order parameter with the 
velocity v which vanishes as v ~ (/ — f c )^ at the transi- 
tion. One also introduces the dynamic exponent z which re- 
lates time and space by t ~ x z . The critical force f c , which 
must be tuned to reach the scale-invariant regime, plays a 
role similar to the critical temperature in thermal phase tran- 
sitions. There are many subtleties however within this anal- 
ogy, since depinning is a non-equilibrium transition at zero 
temperature, where quenched disorder dominates. As the cor- 
responding static problem of elastic manifolds in disordered 
media 1 13], the depinning problem suffers from two peculiar- 
ities, when compared to standard critical phenomena: First an 
infinite set of operators becomes relevant simultaneously be- 



low the internal upper critical dimension d < d uc = 4. Sec- 
ond, their study is more difficult due to the dimensional re- 
duction phenomenon, which renders naive zero-temperature 
perturbation theory trivial, hence useless. The way out in- 
volves first parameterizing the set of relevant operators into a 
function, A(u), the second cumulant of the random pinning 
force. It was shown in Refs. fl4l [l5Tl that the correspond- 
ing functional renormalization group (FRG) provides an ade- 
quate description of the depinning transition if one considers 
the non-analytic renormalized function A(u). It was shown 
only recently that the FRG can be unambiguously extended 
to higher loop order and that the underlying non-analytic field 
theory is renormalizable I la.ll7ll . The FRG equation for A(u) 
has two main non-trivial stable fixed points which describe pe- 
riodic and interface universality classes. Both of them exhibit 
a cusp singularity of the form A* (u) — A* (0) ~ |u| at small it. 
This cusp accounts for the existence of the critical threshold 
force f c ~ A*'(0 + ). The corresponding critical exp onents 
have been computed to second order in e = 4 — d | 16 U17ll . 

Despite this progress many open questions remain. Among 
them is the problem of sample-to-sample fluctuations, i.e. the 
probability distributions of a given observable and their rela- 
tion to extreme value statistics. These were studied mostly for 
static quantities. The distribution of the energy of pinned man- 
ifolds was analyzed in Ref . fl8l fT^l . The distribution of the 
mean squared width of an interface at the depinning transition 
was calculated using a Gaussian approximation for the dis- 
placement field, yielding a result quite close to numerics 1 20] . 
It was shown how systematic corrections can be computed in 
the field theory of depinning within an e = 4 — d expansion 
l2lll . One expects sample-to-sample fluctuations to play an 
important role in the dynamics too, leading to a broad distri- 
bution of time scales. The divergence of the typical energy 
barrier with scale, responsible for the ultraslow creep motion, 
is predicted by phenomenological arguments |7]. A numerical 
study in Ref. 12211 of the distribution of barriers confirms that 
the typical barrier scales as the energy minima, as predicted 
by 1-loop FRG studies l23ll . The more difficult question of 
predicting the distribution of energy barriers was addressed in 
Ref. 12411 using the FRG, and in Ref. [25] using extreme value 
statistics. 

An important and debated question is to characterize the 
finite-size fluctuations of the critical force and whether they 
obey finite-size scaling (FSS). Similar questions were investi- 
gated recently in the context of heteropolymer unbinding tran- 
sitions (the role of critical force being played there by critical 
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temperature) where violation of FSS was found 1 26] . In the 
depinning problem, one difficulty is to define properly the crit- 
ical force and its fluctuations in the limit of large interface (in- 
ternal) size L. A recent and efficient algorithm allows to 
obtain exactly the critical force f c of an interface in a periodic 
medium of period M (i.e. a cylinder), as well as the so-called 
critical configuration. The latter is defined as the last block- 
ing configuration as / is increased up to f c , which also defines 
fc = fc(L, M). One can refer to this definition as an extremal 
configuration in a given sample. The finite-size sample-to- 
sample distribution of f c (L,M) was computed numerically 
112 811 and found to depend on the aspect ratio k = M/L^ of 
the cylinder. This should be expected since for large k one 
recovers a zero-dimensional problem and the interface will be 
blocked by rare disorder configurations, hence dominated by 
extremal statistics. However this results seems to depend on 
the precise definition and one may ask whether a more fun- 
damental definition exists, with no need to specify a value 
for k. An alternative approach is to define the observables 
at the depinning transition as the time average in the mov- 
ing state at fixed velocity v, in the limit v — -> + . This def- 
inition, to which we refer as the quasi-static depinning limit, 
is usually associated to the FRG approach of the depinning 
transition. Observables calculated in this approach must, a 
priori, be distinguished from those computed in the critical 
configuration. Since the time average is usually performed 
in a steady state, to avoid dominance by history dependence, 
it also requires specifying boundary conditions. It is widely 
believed that both approaches give the same result, at least 
for (N = 1 component) interfaces, since in the limit of infi- 
nite systems (L — > oo) all quasi-static configurations should 
have the same statistical properties and the critical force f c 
should be self-averaging. However it is less clear how these 
approaches compare when applied to finite-size fluctuations 
where the dispersion in local pinning forces becomes impor- 
tant. 

In the present paper we study the distribution of the thresh- 
old forces by means of the functional renormalization group. 
Within the field theory we propose two definitions of the criti- 
cal force f c (L) in finite size L and show that they are identical 
to one loop in the renormalized theory. We compute the cu- 
mulants of f c (L) and extract the distribution which is found 
to be universal, up to a shift in / (the critcal force f c ) and one 
scale-parameter (the width of the distribution). All results are 
valid within the e = 4 — d expansion and extrapolations to low 
dimension are discussed. The critical force studied here is de- 
fined from a fixed center of mass ensemble. As we point out 
it can be, in principle, obtained in numerics. Since it does not 
refer to any transverse size M, it is more fundamental than the 
one used in the numerical studies on a cylinder. We discuss 
how the latter one can in principle be recovered. 

The paper is organized as follows. Section UTI introduces 
the model and the FRG treatment of the depinning transition. 
In Section [HI] we compute the bare distribution of threshold 
forces to one-loop order using an improved perturbation the- 
ory, and renormalize it for the case of an elastic interface. In 
Section IIVI we discuss the renormalization for periodic sys- 
tems. In Section [V] we discuss the relation between the dis- 



tribution of critical forces, in the quasi-static limit, and in the 
critical configuration. 



n. MODEL AND FRG DESCRIPTION 

Let us consider the motion of a one-component elastic man- 
ifold with short-range elasticity. The configuration of the 
manifold can be described by a scalar displacement field u x t, 
where x denotes the d dimensional internal coordinate of the 
manifold. We study the over-damped dynamics of a manifold 
in the disordered medium which obeys the following equation 
of motion 

■qd t u xt = cV 2 u xt + F(x, u xt ) + /, (1) 

where rj is the bare friction and c the elasticity. The quenched 
random force F{x, it) can be chosen Gaussian with zero mean 
and variance 

F(x, u)F(x', it') = A(u - u')6 d (x - x 1 ), (2) 

For periodic systems the function A (it) is periodic, while for 
interfaces it decays exponentially for large it. In the latter 
case, in contrast to the statics, at depinning both random bond 
(RB) and random field (RF) microscopic disorder renormalize 
to the same fixed point, which has RF characteristics, so that 
we can restrict ourselves to the latter case, a is the width of the 
function A(u). To make the problem well-defined we imply 
an UV cutoff at scale A -1 . We consider a finite system of size 
L with periodic boundary conditions. The size L serves as an 
IR cutoff i.e. it plays the role of the mass in the corresponding 
field theory. One can easily see that due to the tilt symmetry 
the elastic constant remains uncorrected to all orders so that 
we are free to fix c = 1. 

Below, starting in Section [HI] we will find it convenient to 
work in the comoving frame. To that end we shift u xt — > vt + 

u xt , s.t. (u xt ) = and / -> f-ryo, where v = L~ d {$ x d t u xt ) 
is the velocity of the center of mass. Here the angular brack- 
ets stand for the average over different initial configurations 
(since we are studying zero temperature dynamics) and the 
overline denotes the average over disorder distribution. We 
will assume that a steady state attractor has been reached, 
hence that averages depend only on time differences and not 
on a specific choice of initial conditions. 

To study the dynamics of an elastic manifold efficiently, we 
use the formalism of generating functionals. Introducing the 
response field u xt one can compute the average of the observ- 
able over dynamic trajectories with different initial 
conditions for a particular disorder configuration as follows 

(A[u x<t \) = J V[u]V[u]A[u x , t }e- s ^ . (3) 

Sf is the action for a particular realization of the disorder (a 
particular sample). To compute the average of the observables 
which explicitly depends on the random force at the position 
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of the manifold we introduce the source J x t for the random 
force F so that the corresponding action reads 



Sf[u, u] 



iuxt(rjd t - V 2 )w xt 



iu xt {F(x,u xt ) + f xt } - / J x tF{ )-(4) 

xt J xt 

After averaging over the disorder distribution, any observable 
which depends on the displacement field and the random force 
at the position of the manifold can be computed as follows 



(A[u Xtt \F(x 1 ,u xutl )...F(x n ,u Xnt t n )) 
6 



SJx 1 ,t 1 ---^Jx Tl ,t r , 



V[u]V[u]A[u Xit ]e 



-S[u,u] 



J=0 



S[u, u] is the effective action, which can be split into two 
parts: the free part So being quadratic in fields and the in- 
teraction part S'int containing all non-linear terms 

So= iu xt (r]d t - V 2 )u xt - / iu xt f xt , 

Jxt Jxt 

S'int = -\ ( (iu xt + J xt )A(u xt - U x t')(iUxt> + Jxt')- 
A J xtt' 

Setting Jxt = we recover the action used in Refs. 
The quadratic part So gives the free response 



(u g<t iu- q ,o) = R q .t 



V 



(5) 



while the free correlation function is C q _t = (u q ,t w_ g ,o) = 
at zero temperature. The splitted diagrammatics for the pertur- 
bation theory in disorder A was developed in Refs. lll6lll7ll . It 
is known that naive perturbation theory, obtained by taking for 
A an analytic function exhibits the property of dimensional 
reduction and fails to describe the physics, giving for example 
an incorrect roughness exponent. The physical reason for this 
is the existence of a large number of metastable states. 

Let us briefly sketch the FRG analysis of the system under 
consideration. Power counting shows that the whole function 
A (it) becomes relevant below d uc = 4 and thus one has to 
renormalize the whole function. To extract the scaling be- 
havior one has to study the flow of the renormalized func- 
tion A under changing the IR cutoff towards infinity. Various 
choices for the IR cutoff were discussed in Refs.l 16, d. A 
convenient choice is to add a small mass to, so that the scal- 
ing behavior can be extracted from the effective action T[u, u] 
of the theory as to decreases to zero. To study the finite-size 
distribution of threshold forces, we use, as in [21], the sys- 
tem size L as the natural IR cutoff. Then any integral over 
momentum q has to be replaced by the sum according to the 
rule J — > L~ d Yl q < where the sum runs over all q = 27rfc/L, 

k G Z d , k 7^ 0. Exclusion of the zero mode means that we 
are working in an ensemble of fixed center of mass, a point 
further discussed in SectionlVl 

Let us define the rescaled disorder as 



A(u) = -Li 2C - e A(wir c ) , 
eh 



(6) 



where 1\ = L e I\ = f \q\ 2 is the one-loop integral. It was 
shown in Refs. lll4l Fl51 that the FRG equation, i.e. the flow 
equation for the running disorder correlator can be written to 
one-loop order as 



Ld L A(u) 



(e - 2C)A(u) + C«A'(u) 



(A(«) - A(0) 



(7) 



the two loop flow equation being obtained in Refs. fl6lfT7Tl . 
"0" means a derivative at fixed bare quantities. The flow of 
the correlator is such that A(u) acquires a cusp at the ori- 
gin u = at the Larkin scale L c ~ (c 2 a 2 / A(O)) 1 ^. Be- 
yond the Larkin scale (L > L c ) the renormalized correlator is 
singular and perturbation theory breaks down. Nevertheless, 
the flow tends to a non-trivial fixed-point (FP) solution A* (u) 
with a new value for the roughness exponent which controls 
the large-scale behavior. There are two FPs which describe 
interfaces and periodic systems correspondingly. The former 
FP has C = e/3 + C(e 2 ), while the latter one has ( = due 
to periodicity. The renormalization of the mobility gives the 
value of the dynamic exponent 



z=2+ L—irn 
dL 



= 2- A"(0) + O(A 2 ), 



(8) 



where r/ji is the renormalized mobility. Other critical expo- 
nents can be computed using the scaling relations 



1 



c 



c 



(9) 



To renormalize the theory, one needs an additional counter- 
term for the excess force / — r)v, which comes with an UV 
divergence ~ A 2 . This term is analogous to the critical tem- 
perature shift in the ip 4 theory, and gives the critical threshold 
force /*. It is zero in the bare theory. We expect that in the 
limit of an infinite system L — > oo, the critical force becomes 
sample independent if there is the same distribution of disor- 
der in each sample and thus lirriL—Kxj Pl (/) = 5(f — /*). 
However the situation is different for finite systems. Accord- 
ing to the general theorem for random systems l32ll there ex- 
ists a finite-size scaling correlation length £fs which charac- 
terizes the distribution of the observables in an ensemble of 
samples and which in principle has to be distinguished from 
the intrinsic correlation length £ which enters into correlation 
functions. Approaching the critical point, the finite-size corre- 
lation length diverges similar to the intrinsic correlation length 
as £fs ~ I / ~ fc I ~" FS ■ In general vps is different from v, and 
satisfied the inequality 



^FS > 2/(d + 



(10) 



where d+( is the effective dimension of the disordered system 
considered. Thus for an ensemble of samples of linear size L 
the width of the distribution of critical forces is characterized 
by a scale £fs = L and reads 



{(ML)-f*Y 



L 



-2/Vfs 



(ID 
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FIG. 1: Diagrams Di, (i — 1, 4) contributing to the second cu- 
mulant of the threshold force. We have adopted the splitted diagram- 
matics used in Ref. IrHl : arrowed line indicates a response propaga- 
tor dashed line - the splitted vertex iu x tiu x t' A(u x t — u xt i). The 
corresponding expressions are given by Eq. 1 19i . 



For periodic systems £ = and v = 1/2 so that vps ^ v 
for d < 4. For interfaces it was proposed 11511 that vps = v - 
While z^fs satisfies d!0> with the equal-sign in 1-loop order, 
at 2-loop order the inequality becomes strict. As discussed in 
Refs. fl5l[l7ll this difference is closely related to the stability 
of the FP which controls the scaling behavior. 



III. DISTRIBUTION OF THRESHOLD FORCES 

A. Perturbation theory 

Let us show how the critical force distribution can be com- 
puted within "improved" perturbation theory. "Improved" 
means that we assume the the disorder correlator A(u) to be 
non-analytic, since for analytic disorder perturbation theory 
gives a zero-threshold force. Then using FRG we renormalize 
our result to 1-loop order. Effectively, this is a summation of 
an infinite subset of diagrams. We define the distribution of 
threshold force as follows 



PL{fc) = (s(f c + L-*J F 



(x,vt + u xt ) 



(12) 



From now on we work in the comoving frame and the average 
is performed with the action S in the quasi-static limit v — > 
+ , as is usually done in the FRG theory of the depinning 
transition. 

Let us introduce the corresponding characteristic function 

P £ (A) - (e-^W) = J df c e~^P L (fc) (13) 



which can be expressed through the cumulants (/ c (L) n ) as 
follows 



p L (A) = ex P hr 



HA)" 



(fc(LY 



(14) 



The computation of the first cumulant to one loop is trivial. 
The random force that the interface actually feels in the point 
x is given by 



(F(x,u x>t + vt)) = ( / A(u xt - u xtl + v(t - ti))iu xtl 



A'(v(t - t 1 ))[R a -o, t -t 1 - Px=o,t=o] 



(15) 



We will adopt Ito's prescription in which R x .t=o = 0. Note 
that this corresponds to the definition O(0) = 0. Taking the 
quasi-static limit v — ► + , we obtain the well-known expres- 
sion for the average critical force 



ft = - A '(o + ) / 



dtR T 



A'(0+) 



(16) 



Note that the critical force (I16t diverges at large momentum 
as A d ~ 2 and therefore is not universal, i.e. it depends on mi- 
croscopic parameters. This is analogous to to the shift of the 
critical temperature in standard critical phenomena, caused by 
fluctuations. As we know this is shift is also non-universal. 
However, we expect that the distribution of critical forces for 
a finite system around the average value are universal, once 
the distribution is properly normalized. The computation of 
the the n-th cumulant is more tricky. Before considering the 
general case let us show how this works for the second cumu- 
lant. Using the generating functional, we can write down the 
formal expression for the effective force-force correlator 



(F(xi,u Xut + vt)F(x 2 ,u X2i t + vt)) = A(0)5 d {x 1 - x 2 ) 



A(u xit 



A(u X2 



+ v(t - tx))iu xitl 



2t 2 + v(t ~ t 2 ))iu X2ti 



(17) 



The first term on the r.h.s. of Eq. dl7l is the bare disorder dis- 
tribution. It is given by Eq. and is a pure Gaussian distribu- 
tion with zero mean. However the moving manifold explores 
a different distribution that is an effective distribution which 
one can observe "sitting" on the moving interface. The second 
term on the r.h.s. of Eq. d!7t as well as the mean value dl 5t are 
the deviation of the effective distribution from the bare one. 
Only connected diagrams contribute to the second cumulant. 
Integrating the second term in Eq. (I17> over fields with the 
weight e _s we obtain the four connected diagrams shown in 
Figure [J The corresponding expressions can be rewritten as 
follows 



tl«2 



A'(v(t — ti))[R X2 - Xl ^-t 1 — Rx 2 -x 1 ,t 2 -t 1 ] 

A'(v(t — t2))[R Xl - X2 ,t-t 2 ~ Rx 1 -x 2 ,t 1 -t 2 



(18) 



To compute the contribution to the variance of the critical 
force, we have to integrate over x\ and x 2 and then multi- 
ply by L~ 2d . This computation is more convenient in Fourier 
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representation: 

D 1 + D 2 + D 3 + D 4 = L- d [A'(0+)] 2 



* I [Rq,t-ti 
'8*1*2 



R, 



g,*2— * 



l][Rq, 



q,t—ti 



R, 



q,ti—t2 



I- (19) 



Due to causality we have D4 — 0. The other diagrams read 

- d f [A'(0+)] 2 



Di = -D 2 = -D 3 = L- 



Summing all contributions we obtain 



(<7 2 ) 2 



(fc(Ly 



L- d A(0)-L- d 



[A / (0 + )] s 

2\2 



(<z 2 ) 



(20) 



(21) 



(/ c (i)») c = (-lT(n-l)\L- d ^ /[A'(0+)] 



We have also derived Eq. i2l\ by using direct perturbation 
theory instead of the generating functional. 



The above calculation can be generalized to arbitrary n. It 
can be simplified significantly if one takes into account that 
all intermediate times ti must be smaller than the observation 
time t: ti < t {i — 1, ..,n). For the n-th cumulant (n > 2) we 
have 



/ [Rq,t-t! — Rq,t 2 -ti] [Rq,t-t 2 - Rq,t 

Jti.. .*„ 



[Rq,t-t n _ 1 — -Rg,t„-*„-i] [R 



q,t-t„ 



R q , tl -tJ- (22) 



Here the factor (n — 1)1 results from different contractions of 
u Xi t - u Xiti and u Xjtj = 1, . . . , n) that form a closed 
loop. Expanding the integrand in Eq. \22\ we find that all 
terms gives the same contribution up to a factor of ±1, except 
for the term composed only from the second response function 
in each bracket. This term gives a closed loop of response- 
functions, which is zero by causality. Using the identity 



where we have taken into account that A'(0 + ) < 0. 

As follows from the above computation, the distribution of 
the critical force can be related to the effective action T[u, u] 
which is a generating functional for one-particle irreducible 
(1PI) vertex functions rV^;^ with S external fields u and 
E external fields u 



(23) 



i=0 



where C % n is a binomial coefficient, we can simplify Eq. (12 2 1 
to 



-(n-l)\L- d ^ I [A,( "| )r . (24) 



We are now in a position to construct the characteristic func- 
tion 



lnP(A) 



1 



A{0)L- d X 2 



9n=l 



(-1)™ /A'(0+) 



L- d i\ 



q- 



(25) 



The latter is nothing but the Taylor series of the logarithm, 
which allows to rewrite Eq. d25i as 



P(A) 



exp 



1 



d \2 



A(0)L- d X 



(26) 



n^n 



Sua 



-T\u, ii] 



(27) 



u—u—0 



Indeed as already seen from the bare action the average thresh- 
old force can be expressed as vertex function r^ 1,0 ^ = 
0,ui = 0) in the quasi-static limit v — > + . Analogously the 
higher-order cumulants can be identified as the higher-order 
vertex-functions according to 



(ULY) c = L-(»-i>* rfc2({ ft = 0, Ui = 0}). (28) 



The general properties of vertex functions d28i for even n was 
discussed in Ref. El . In particular it was noted that loop 
diagrams that contribute to the vertex r' 2 "' ' precisely can- 
cel each over so that the result is given by minus the missing 
contribution from acausal loops. It is easy to verify by direct 
inspection of the Feynman diagrams that definitions d!2t and 
d28t i give the the same result in the one-loop approximation, 
but the question of their equivalence to all orders is open. 
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FIG. 2: Example of diagrams contributing to the third cumulant of 
the critical force. To renormalize these one-loop diagrams to the low- 
est order we replace the disorder lines by vertices T^r (highlighted 
on the left diagram) and (highlighted on the right diagram) 

which are defined in Eqs. 1351 and 1361 and depicted in details in 
Figure [3] 



1 a t 



* t 



r 



(+) 



q - 



1 * h 



i h 



r 



(-) 



FIG. 3: Vertex function F^J (t, ti,t2\q,—q) at the tree level. The 
functions are distinguished by whether the line entering at ti and 



outgoing at ti is "closed" TV^, or "open" T^t^. In the first case 
the time-ordering along the loop is continuous: t\ <ti, while in the 
second case it is interrupted, i.e. there are no restrictions on t\ and ti 
(see also Figure[2j. 
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B. Renormalization 

In this section we focus on the interface problem, periodic 
systems being considered in the next section. The distribu- 
tion of the critical forces in (I26> has been obtained from the 
"improved" perturbation theory, and thus, it cannot be repro- 
duced within the Larkin type models in which all observables 
depend only on A(0). However in the bare theory the disorder 
correlator is an analytic function so that to make the calcula- 
tion consistent we have to first renormalize our theory. To that 
end we replace the bare correlator by the renormalized one. 
This can be viewed as a partial summation of an infinite series 
of diagrams. If we want the distribution of f c strictly to lowest 
order in e = 4 — d then the work is essentially done. However, 
we will demand a bit more and take an additional effect into 
account: replacing the bare correlator by the running one in a 
particular diagram we have to be careful because the scale de- 
pendence acquired by the correlator may be in form of either 
dependence on a mass (as in quantities which do not contain 
integration over momentum) or dependence on the loop mo- 
mentum (which has to be integrated out), or a combination 
of both. The final result presented here will thus be exact to 
lowest order in e and in addition will contain some effects be- 
yond that order (although a full fledge two loop calculation is 
not attempted here). This will allow us to discuss in the next 
section some extrapolations to low dimension. 

Let us start from the renormalization of the first cumulant, 
i.e. the average critical force d!6l >. We remind the reader that 



the average critical force is a non-universal quantity and af- 
terwards we will subtract it and consider the shifted distri- 
bution which is a universal function. According to Eq. l|6} 
the renormalized disorder correlator acquires in the vicinity 
of the fixed-point a scale dependence. Integration over scales 
beyond the Larkin scale yields (see Ref |23, 39] for details): 



</ c (L)) c «-A*'(0+) 



A2-C 
2-C 



(29) 



where in this formula the UV cutoff A is meant to be the 
minimal length of pinned segments of the manifold, i.e. , the 
Larkin length A ~ L~ . 

We now consider the renormalization of the second cumu- 
lant (variance). The corresponding bare expression can be ex- 
pressed through the 2-point vertex function as follows 



A(0) - f 

J a 



A'(0+) 2 



0(A 3 ) 



= L- d r£% = 0;g = 0). 



(30) 



As it was shown in Ref. Jl7l I2TII - the 2-point vertex function 
does not depend on times and scales, 



r%(q) = L^J-A*(0)F 2 ( q L), 
he 



(31) 



with F 2 {z) = Bz £ ~ 2C - + 0(\nz/z 2 ) for large z and F 2 (0) = 
1. Note the the constant B depends on the IR cutoff scheme 
l2lll . Combining Eqs. POt and OH we obtain 



he 



(32) 



We note that this is consistent with the finite-size scaling pre- 
diction: 



(33) 



using v = 1/ (2 — . As we will show below, the full (shifted) 
distribution is also consistent with this scaling. 

To proceed further, let us first consider some typical dia- 
grams contributing to the third cumulant of the critical force, 
which are shown in Figure [2] To renormalize them at lowest 
order, we have replaced the disorder lines by the three-point 
vertices defined as follows 

r ifi« (*>*!> *a; 9i>«2) = r^(t,*i s t a ;«i,ft) 

+ ri^(Mi,t 2 ;<?i,<&). (34) 

At tree level, the vertex function (13 4> can be expressed by di- 
agrams shown in Figure [3] and the corresponding expressions 
read 

r ( Mu(t,tl,t2;qi,q2) = A'(0+)sign(t-t 1 )^-t 2 ),(35) 
^ u (t,ti,t 2 ;qi,q 2 ) = A'(0 + )sign(t 1 -i)<5(i 1 -t 2 )(36) 

Then the summation of diagrams contributing to the the n-th 
cumulant with n > 2 can be carried out along the lines used 
for the bare cumulant and gives 
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</ c (L)«) c = (-1)>- l)\L- d ^ [ [ [ \r^(t,t 2 ,n;q-q)R q , Tl - tl 

r ifi«(M3,'T2;<7,-<?)-R g ,T2-t2 ■■■ r i 2 ^ ) (i^n ) T «-i;'?:-'7)-R<z,r„- 1 -t„- 1 (t , t x , t„ ; q, -q)R q , T 



(37) 



Substituting the tree-level expressions J34I-J36I to Eq. ( I37t 
we recover the bare cumulant i22\ . 

In Appendix 1X1 we compute the vertex function T*t) to 
one-loop order and obtain its large-q asymptotics which reads 



(38) 



where times t\ and t are taken infinitely apart (hence it is a 
quasi-static quantity). The amplitude A and the exponent ip 
are given by 



A = 1 A'*(0+)(1 + O(e)), 
eh 



(39) 



(40) 



and we argue that rp is a new exponent (see Appendix 1X1. 
Taking into account this momentum dependence in Eq. (I37> 
should result in an improved renormalization scheme (com- 
pared to simply replacing the bare quantities by g-independent 
but scale-dependent renormalized parameters) with a different 
form for the q summations appearing below. 

After substituting Eq. (1381 in Eq. d34K the integration over 
times can be performed in the same way as for the bare cumu- 
lant and gives for n > 2 1 40 ] : 



</ c (L) n > c = -{n-l)\A n L n ^ 
xL d 1 ' 



, {qL) n V-^ ' 



(41) 



Note that A ~ A*'(0 + ) < 0. To construct the characteristic 
function let us redefine A -> (2n) 2 ^ X/ \\A\L^ 2 ), that cor- 
responds to measuring f c in units of |^4|L^~ 2 / (2tt) 2 ~^. This 
is a non-universal scale as the value of A*'(0 + ) is not uni- 
versal at the depinning transition. However as we now show, 
once rescaled the (shifted) distribution is universal. 
The characteristic function can be written as 



2" " ^ \ |k| 2 "^ 

keZ d ,k#0«=3 



(42) 



where 



(2tt) 4 A*(0) 
A 2 he 
A*(0) 



f/lTT 



[A'*(0)] 2 



(2ir) 4 - 2 ^ = — (l + 0{e)). (43) 



Gtt 



To compute this universal ratio to lowest order in e we have 
used the 1 -loop FRG fixed-point equation evaluated at u = 0, 



i.e. (e-2C)A*(Q) = A*'(0+) 2 and used ( = e/3 + 0(e 2 ), a 
well as eli — l/(87r 2 ) + 0(e). Summing over n we obtain 



lnP(A) 



\a 2 X 2 



E 

keZ d ,k^0 



iX- 



1 



2 |k| 2 ( 2 -^) 



In 1 



|k|2-V- 



(44) 



C. Lowest order in e = 4 — d 

To obtain the distribution within the e expansion it is suffi- 
cient to set i/j = in the formula above, and to compute all 
sums in d = 4. Let us first give the skewness and kurtosis to 
lowest order in the e expansion. One uses l2lll : 



y - 

£—1 If : 



1 



keZ d ,k^0 



k| 2 P (p-1)! 



y — 



/ dtt p - 1 [6(3,0,e- t ) d - 1], (45) 
Jo 



keZ 4 ,k^0 



kl 6 



E - 

^ If s 



14.8298, 



10.2454, 



(46) 
(47) 



keZ 4 ,k#0 



where 6(3, 0, e *) = X^fcez e is the elliptic theta func- 
tion. Hence: 



03 



0-4 



(/-/) 3 2 



E - 



k£Z d ,k/0 



kl 6 



(48) 



0.0650861e 3/2 , (49) 

^4 „4 2^ IH8 W 



a- — |k|'' 

keZ 4 ,k^0 



-0.01753e^ 



(51) 



Next one can resum to obtain the characteristic function to 
lowest order in e: 



lnP(A) = -^<r 2 X 2 - F(-iX), 



(52) 



F(-iA) 



°°y(e 4At -l-zAt+iA 2 i 2 ) 

x[9(3,0,e-*) 4 - 1]. (53) 



This result can be reexpressed as follows. As d — > 4 
the shifted dimensionless critical force / = (f c (L) - 
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f c (L))/-J fc(L) 2 becomes a univariate gaussian random 
variable. For small e > 0, f can be (formally) expressed as the 
sum of two independent random variables / = /o + 
where /o is gaussian of variance 1 + 0(e) and /i is a ran- 
dom variable of order unity with a non trivial distribution, the 
logarithm of its characteristic function being given (up to a 
quadratic term) by F(— iX). 

We now analyze the shape of these distributions in physical 
dimension. 



D. Fourier inversion 



-15 



-10 



0/08 



0.06 



0.04 



0.02 



/ 



PAf) 



10 



15 



In this section we compute the inverse Fourier transform of 
Eq. (1441 in physical dimensions, using our improved scheme. 

Let us start by discussing d = 1. We use a natural extrapo- 
lation of our above result, setting e = 3 in the above formulae 
(which are exact to lowest order in e). From Eq. J44i we ob- 
tain 



FIG. 4: The three indistinguishable curves for shifted critical-force 
distribution (d = 1): improved perturbation theory (ip = 0), renor- 
malized to one-loop theory and approximate expression I56i . Here 
a — 7ry2. The difference between distributions is shown in Fig- 
ure|5l 



P(A) = cxp[-A 2 ]n 1- 



fc=i 



iX 

P/3 



x exp 



iX 



X - 



1 



£2/3 £.4/3 



(54) 




where we have used -0 = 4/3 and a = n y/2. The inverse 
Fourier transform computed numerically is shown in Figure^] 
Eq. <!44i with ip = can formally be considered as the result 
of improved perturbation theory in non-analytic disorder. The 
inverse Fourier transform of the latter is also shown in Fig- 
ure |4] and can not be visually distinguished from the shifted 
renormalized distribution. The renormalized distribution is 
more appealing since it guarantees a non-negative defined crit- 
ical force, which is not the case for the bare one, especially in 
d = 1 where the bare averaged critical force dl6l is finite. By 
contrast the averaged renormalized critical force is controlled 

by the UV cutoff and is of order L~ 1,v while the typical flue- critical force distribution (04) can be written as 
tuation is much smaller, of order L -1 /", in the limit of in- 
terest, here L ^> L c . It is interesting that the renormalized 
distribution is well approximated by the Eq. J54i in which we PW = ex P 

keep only the first factor with k = 1: 



FIG. 5: (Color online) Shifted critical-force distribution for d — 1 
with Gaussian subtracted. The dashed lines: approximate expression 
1561 and improved perturbation theory (tp = 0); solid line: renor- 
malized to one-loop theory. 



1 ^2 \2 

2 



P(X) w exp 



\2 

2 



{(1 - iX) 2 exp [2iX - A 2 ] j , (55) 



x exp 



iX- 



n 

k£Z d ,k^0 



1 - 



iX 



k| 2 -</> 2 |k| 2 ( 2 ^) 



(57) 



The inverse Fourier transform of Eq. d55l > reads 



Analogously to the case d = 1 the inverse Fourier transform 
of Eq. J57> can be well approximated by 



P[f] 



((^ + / + 4) s 



(2 + ( j 2 ) 5/2 v / 2^ 



(2 + /) 2 
-g 2(2 + ^^) 



(56) P(X) = exp 



-a 2 A 2 



{ (1 - iX) 2d exp [2diX - dX 2 ] } (58) 



and is also shown in Figure |4] The difference in the critical 
force distributions obtained within improved perturbation the- 
ory, renormalized to one-loop theory and approximation d56t 
is indicated in Figure [5] 

For a d-dimensional system the Fourier transform of the 



which does not depend on if;. 

We now compute the standard deviation, and the kurtosis of 
the above distributions. The standard deviation reads 



(/-/) 2 =^. 



(59) 
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The skewness is defined as 
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U-ff 



^3 |k|3(2-V) 



(60) 



For d = 1 we obtain a 3 = V2/(6tt) ~ 0.075 [0.046 for the 
distribution 1561 1. The positive value for the skewness indi- 
cates that the right tail of the distribution is heavier than the 
left tail. 

The kurtosis for a Gaussian distribution is three. For this 
reason, excess kurtosis is defined as 



(/ - /) 4 „ 3! 

CT 4 = 1 i = 7 



E 



1 



| k |4(2-V) 



(61) 



For d = 1 we obtain a 4 = -3C(8/3)/tt 4 ~ -0.04 [-0.031 
for the distribution ( I56H . Here ((x) is the Riemann Zeta 
function. The small negative excess kurtosis indicates that 
the distribution is slightly more flat than a Gaussian distri- 
bution, while the deviation from the Gaussian distribution is 
quite small. 



IV. PERIODIC SYSTEMS 

The renormalization of the critical force distribution for pe- 
riodic systems requires a separate consideration. Indeed as 
was shown in Ref. 11511 . in the periodic case there is an ad- 
ditional relevant operator which is the uniform part of A(u) 
so that the RP fixed point is unstable. The flow equation for 
this operator can be derived by integration of the RG equation 
over one period fl7ll 

l l 
Ld L J A(u)du = e J A(u)du + 0(A 3 ), (62) 
o o 

where we have explicitly used ( = 0. Thus in the vicinity of 
the RP FR the flow of the dimensionless disorder is given by 

A(u) = A*(u) + cL £ (63) 

where the non-universal constant c can be estimated as fl7ll 
l 

c = L~ e J (A( barc )(u) - &*(«)) du 



-L-z I A*(u)du = I-* ( — + 0{e) ) > 0. 



This runaway correction to the scaling behavior at the RP FP 
contributes to all quantities which depend on A(0) but not 
to those which depend on A'(0). Therefore in the case of 
a periodic system the renormalized second cumulant of the 
critical force reads 

1 



(/ c (i) a ) c = ^A*(0)L- 4 + -^L- d . (64) 
he he 



Higher order cumulants are still given by Eq. (14-11 with Q = 
and therefore scale with L as (f c (L) n ) = L~ 2n . We would 
like to emphasize that the correction to scaling in Eq. (I64> de- 
scribes the sampel-to-sample fluctuations and can not be seen 
in one sample because in each sample there is only one pinned 
configuration. As a result for d < 4 only the second cumulant 
of the scaled critical force, which scales as L~ d , survives in 
the limit L — > oo resulting in vps = 2/d and a pure Gaussian 
distribution for the scaled critical force. 



V. DISCUSSION 

In the present paper we have computed the renormalized 
distribution P[f — rjv] averaged over all pinned configura- 
tions in the limit v — > + , which we identify with the crit- 
ical force distribution Pl(/ c ). The average critical force is 
a non-universal quantity which depends on microscopic de- 
tails of the interactions like the UV cutoff as well as on de- 
tails of the disorder distribution. After subtraction of the aver- 
age value, the shifted distribution of the critical force contains 
only one non-universal scale which can be fixed, e.g. by fixing 
the second cumulant. The resulting dimensionless distribution 
is then fully universal, i.e. it does not depend on properties at 
small scales. We have computed it taking into account only 
the second cumulant of the bare disorder distribution, since 
it is the only cumulant relevant in the RG sense. Higher cu- 
mulants, which are generated by coarse graining are irrelevant 
operators and their contribution to the cumulants of the critical 
force must carry additional dependence on the cutoff. Hence 
we expect that they result only in a shift of the non-universal 
expectation values. 

Let us now discuss the role of the transverse sample size M 
(size of the box). In numerical studies of depinning of elas- 
tic interfaces, either via an exact determination of the critical 
state, or via Langevin dynamics 112(1 1271 l28l l29t 13011 a cylin- 
drical system which is periodic in both directions was stud- 
ied: longitudinal with period L and transverse with period M. 
Hence this is equivalent to a periodic disorder with period M. 
It is known that a periodic system has a unique pinned config- 
uration for any period M [38]. If M <C L, this configuration 
spreads out through the whole box L d x M and there is only 
one independent pinned configurations. As we have shown for 
the RP class the distribution of critical forces is Gaussian, and 
thus, for elastic interfaces in the limit L — > oo, with M fixed, 
the distribution of the critical force also becomes Gaussian. 

The case where the period M is taken to infinity at the same 
time as L is relevant for elastic interfaces and quite different. 
The pinned interface has a r.m.s. width w — k w L^ so that in a 
sample of transverse size w it also has one unique statistically 
independent pinned configuration. One may then argue that its 
critical force distribution is Pl(/ c ) whose characteristic func- 
tion is given by Eq. J44I . In the numerical studies one should 
thus be careful in choosing the size of the periodic box M. 
If M scales like L 1 * with (' < £ the system will crossover 
from the random manifold to the RP FP and the finite-size 
scaling analysis will results in some mixture of interface and 
periodic system properties, while the critical force distribution 
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will tend to a Gaussian one. On the over hand if M 3> w, the 
sample can be divided into about M/w subsamples, which 
can be argued to be (almost) statistically independent, with 
independent pinned configurations. Each configuration has 
a slightly different critical force which is distributed accord- 
ing to our FRG result. If one defines the total critical force 
as a maximum of all the critical forces of these subsamples, 
it becomes A/-dependent and its shifted distribution tends to 
the distribution of the extreme value statistics [28]. Following 
Ref. 12811 let us introduce for every configuration a of the in- 
terface in the sample of width M = fci> the depinning force 
fd(oc) and then associate the threshold force of the whole sam- 
ple with the following maximal value /£ = max a {fd(a)}. 
In each sample there are only w M/w independent pinned 
configuration, so that the distribution of the maximum of the 
corresponding critical forces can be written as 



p M (r ciM /w) = — 



d 



fc 



df'P L [f] 



M/w 



(65) 



According to the general theorem of extreme value statistics 
l33ll . for large samples, i.e. , in the limit M/w — > oo this 
distribution approaches the Gumbel distribution. The latter is 
provided by the tails of the distribution of the critical force for 
each independent pinning configuration, as given in Eq. {SJ. 
According to this distribution the average maximal threshold 
force of samples of size M behaves as \n(k/k w ). For large 
samples with M ^> w it can be extremely large. The above 
procedure completely washes out all details of the underlying 
distribution Pl(/) computed here, except for its width, and 
replaces it by the model-independent function obtained from 
extreme value statistics. As an illustration, we have plotted 
the force distribution obtained using d65l > for M/w = 10 on 
figurelA5l 

The above arguments suggest that the critical force distri- 
bution computed here via the FRG should be compared with 
the numerics on a cylinder of aspect ratio parameter k sa k w 
defined above from the width. We now propose a more pre- 
cise statement to identify the critical force computed in this 
paper. We note that in the calculations performed here within 
the FRG the position of the center of mass was held fixed 
(since all momentum integrations excluded the uniform mode 
q = 0). Hence we are working in the fixed center of mass 
ensemble. This suggests the definition: 

f c {u ,L) = max {f d (a(u Q , L))} (66) 

a(u .L) 

where the maximum is over all configurations a(uo, L) with 
center of mass uo and length L (and periodic boundary condi- 
tions along the interface). It can in principle be evaluated nu- 
merically by direct enumeration for a discrete interface model. 
One can then check that it has a well-defined L — > oo limit 
with no need for a transverse box, and one can then numeri- 
cally compute the finite-size distribution for the ensemble of 
a(iiQ, L). This distribution should identify with the one com- 
puted here within the FRG (in the massless scheme). It is a 
more fundamental object than the critical force defined on a 




/-/, 



FIG. 6: (Color online) The normalized critical-force distribution for 
d = 1. The solid line is the distribution for the interface in the box of 
size w given by Eq. I56i . The dashed line is the Gumbel distribution, 
i.e. the distribution of the maximal threshold force in the limit of an 
infinite box. The points are computed using Eq. 1651 for the interface 
in the finite box of size M — Ww. 



cylinder. The latter can then be retrieved in principle as 



max/ c (u ,L) 

"0 



(67) 



on the same cylinder, leading to extremal statistics, as dis- 
cussed above. 

The above considerations illustrate that the statistics of the 
depinning threshold force at finite-size is a rather subtle ques- 
tion. Many questions remain open. It would be interesting to 
find the proper steady state corresponding to the above defi- 
nition ( I66> . Also a systematic study of memory effects in the 
threshold force would be of high interest especially regarding 
experiments. Indeed, these memory effects may be of impor- 
tance for aging 1 36] and hysteresis phenomena 1 37] controlled 
in some regimes by the slow dynamics of domain walls. There 
the observed threshold force may not be the largest one but a 
threshold force which characterizes a piece of the sample in 
which the interface got trapped. Hence the data should be in- 
terpreted with care to disentangle history effects from finite 
size effects. It would be very interesting to develop numerical 
schemes to investigate these questions in particular an efficient 
algorithm to compute J66l >. 
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FIG. 7: 1-loop dynamical diagrams correcting A. 



APPENDIX A: CALCULATION OF THE ASYMPTOTICS OF 

r (+) 

uuu 



To renormalize the n-th cumulant of the critical force we 

-.(+) 



and 1\ — I\ (0), the contributions to JA1I are as follows 

- { [h (q)] + [h (0) + h (q)\ + [0] + [h (0) - /! (0)] } 

xA^(O)A^O) 
+ {[-h(0)] + [0] + [Jx(0)] + [0]}A (0), A-(0) (A4) 

where terms in rectangular brackets are in the order of their 
appearance from diagrams a, b, c and d of figure [7] We re- 
mark that contributions proportional to Ao(0)Aq"(0) exactly 
cancel, and we obtain dAli . 

To renormalize the vertex function flAlb we have to reex- 
press the bare disorder correlator by the renormalized dimen- 
sionless one: 



A (u) = 
m e (a(u) + [A'( 



(A(ti)-A(0))A"(tt)]/i}. 



need the leading asymptotics of the vertex functions TVt^ and Differentiating the latter expression w.r.t. u we get after 



T^J. The latter is insensitive to the IR cutoff scheme up to 
corrections to the amplitude. In order to extract this asymp- 
totics, we adopt a massive scheme for the IR regularization 
because it leads to the simplest calculations. If one then needs 
the corresponding amplitude in a different scheme, one can 
relate it to the amplitude in the massive scheme following the 
methods developed in Ref . H35I1 . 

Straightforward perturbation theory gives to first order in 
the bare disorder denoted here A for the "open" vertex TVV 



Tit ) u (t,t 1 ,t 2 ;q,-q) = A' o (0 + ){l-A'^0) 



+ 



(p 



(p 2 + m 2 )[(p 
1 O 



-,2\2 



q) 2 + m' 

-p 2 (t-ti) 



(Al) 



-q) (A2) 



/ rU(Vi,* 2 ;<?,-?) = - / ^ u (t,ti,h;q,-q) 
Jt 2 Jt 2 

Identity (IA2> holds to all orders by definition. The last 
term in Eq. JA1 1 reflects the dynamic nature of the vertex 
rvtL(i 3 t\, t2\ q, —q). However, if we integrated this term also 
over ti the result will not depend on the observation time t. As 
a consequence, the n-th cumulant of the critical force is deter- 
mined by the quasi static integrals (to any order in disorder). 

The easiest way to obtain (IA1I is as follows: One starts 
from the 1-loop diagrams given in Figure [7J expands to first 
order in the fields u, and then replaces Aq™^ (itt — Uf) by 
Aq (0 + ) x [sign(t — t')} n . The diagrams are then split into the 

two classes and T^t^, depending on whether the single 
field u is connected to u(t\) or ufo). In order to extract the 
time-independent terms, one choses t — t% — > oo. This pre- 
scriptions allows for an easy integration of all 16 diagrams. 
Denoting 



h(q) 



1 



p (p 2 + m 2 )[(p + q) 2 + m 2 



m E Ii(q/m) , 
(A3) 



rescaling A(u) = -i-m 2 ^A(™^) 



A o (0) 



^A'(0+) 
eh 

1 



1 



l + ^-A"(0)m ! . 

eh J P {P 2 + m 2 ) 2 



A^'(0) = -^m £ A"(0 + )+O(A"(0) 2 ,A'(0 + )A'"(0+) . 

eli \ ' 

Omitting the last term in Eq. JA1I we obtain the following 
expression for the renormalized vertex function 

riti(Mi,i 2 ;g,-g) = ™ £ - c ^A'(0+) 

eh 

\-2h"($>)\rrf[h{q)-hm 
eh 

Note that depinning of the non-periodic systems is described 
by the fixed point with A 
integral h (y) reads ifTvtl 

l 

da 



by the fixed point with A"(0)* = |e + 0(e 2 ). The one-loop 



(A5) 





where K d = 2ir d / 2 / ' ((2i:) d T {d / 2)) is area of a d dimen- 
sional sphere divided by (27r) rf . Taking into account that 
h = h(0) = J q {q 2 + I)- 2 = \K d T (|) r (f) we obtain 

ri+l(Mi,t 2 ;<z,-<z) = m^ c 4 A '*(° + ) 

eh 



da In [1 + a(l - a)y 2 ] + 0{e 2 ) 



We are interested in the asymptotic behavior for z — * oo. In 
this limit we have 



da In [l + a(l - a)y 2 ] = -2 - 



In 2 



- ln(2 + y 2 - y^A + y 2 )} = -2 + 21ny + ( ^ 
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Matching to a power-law asymptotic behavior we find Replacing m by 1/Lwe obtain Eq. J38I . 

l + ^e(-2 + 2lny) + 0(e 2 )^y^ 9 (l - ~ £ ) . (A6) 
As a result we obtain for g/m ^> 1: 

/ r«(t,*i,* 2 ; <?,-<?) 
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